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Abstract 

We consider the superposition of the first two members of the gravitational hierarchy (Ein- 
stein plus first Gauss-Bonnet(GB)) interacting with the superposition of the first two members 
of the SOr±\ (d) Yang-Mills hierarchy, in d dimensions. Such systems can occur in the low energy 
effective action of string theory. Particle-like solutions in dimensions d = 6, 8 are constructed 
respectively. Our results reveal qualitatively new properties featuring double-valued solutions 
with critical behaviour. In this preliminary study, we have restricted ourselves to one-node 
solutions. 
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1 Introduction 



Gravitational theories in higher dimensions are of current interest in the contexts of the AdS/CFT 
correspondence and of theories with large and infinite extra dimensions with non-factorisable met- 
rics. Since non-Abelian gauge fields feature in the low energy effective action of string theory, it 
is interesting to study the properties of the corresponding Einstein- Yang-Mills (EYM) systems. It 
is our purpose in the present work, to study higher dimensional EYM systems extended by higher 
order terms in both gravitational and YM curvatures. 

In particular, we construct static particle-like solutions in the higher curvature version of Einstein- 
Yang-Mills systems in ci-dimensional spacetimes in which the SO(d) gauge field takes its values in the 
chiral representation SO(d)±. This is analogous to the 4 dimensional model with SO ± (4) = SU±(2) 
Yang-Mills (YM) field interacting with Einstein gravity, in the work of Bartnik and McKinnon [|l| . 

In the original system studied in the gravitational action provides the required scaling be- 
haviour needed to balance the scaling of the static Yang-Mills (YM) system in 3 Euclidean dimensions. 
This scaling argument is not rigorous because the gravitational action is not positive definite by con- 
struction, but nonetheless it enables the construction of particle like solutions. As a result these 
solutions are unstable and can be viewed as sphalerons|2|] rather than solitons. The (d — 1) static 
solutions we construct here are expected to have such an instability, though in this preliminary work 
we will not carry out the corresponding stability analysis. We will restrict to the simplest model 
incorporating only up to the second order curvature terms in both the YM and gravitational sectors. 

As it happens, to go to 8 > d > 4 it is necessary to include the second member of the YM 
hierarchy [RJ to provide the requisite scaling, even in the absence of the gravitational GB term. 
(Similarly for 12 > d > 8 it is necessary to include the thrid member of the YM hierarchy.) We 
will restrict our considerations to the first two members of the YM hierarchy, and hence to d < 8. 
Thus we will consider systems (a) featuring only Einstein gravity, (b) only GB gravity, and (c) mixed 
Einstein and GB gravity, but always with the first two members of the YM hierarchy. 

In section 0, we introduce the models to be studied and in section | we derive the classical 
equations subject to our spherically symmetric Ansatz. This consists of the rf-dimensional analogue 
of the 4-dimensional metric Ansatz employed in |l|, , and the spherically symmetric Ansatz for the 
YM field in the SO±(d) (chiral) representation in <i-dimensions, analogous with the Ansatz for the 
SU{2) YM field in the SO±(4) (chiral) representation in 4-dimensions used in [J. Section f| and 
its subsections are devoted to our numerical study of these systems and we give a summary and 
discuss our results in section |^. 



2 The model(s) 

The generic system of interest in this work is the superposed Yang-Mills hierarchy || interacting 
with the superposed gravitational hierarchy consisting of the Einstein and all possible Gauss-Bonnet 
(GB) terms in a given spacetime dimension d. 

£ = £grav + £-Y M ■ (1) 

In this preliminary work, we will carry out this construction for dimensions d = 6, 8 only. 
The definition we use for the superposed gravitational hierarchy is 

Pl 1 

£grav = X/ o"" K P e ^(p) ' (^) 
p=l *P 
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in which the p = 1 term describes the usual Einstein gravity, and the higher p terms the successive 
GB contributions. In (Q) e = det = ^Jdet g^, are the vielbeins, and is the p-Ricci scalar. 

In terms of the totally antisymmetrised p-fold product of the Riemann tensor Wz, yielding the 2p- 
form p-Riemann tensor (with the notation [a&c...] implying total antisymmetrisation of the indices 
a, b, c, ...) 

T}aia,2...a2 p _ E>[ a l a 2 00304 D i2p-ia2p] /o\ 
^Vllta— A*2p \ptlH2 M3At4'--- n 'M2p-lA12p] ' V / 

we define the p-Ricci scalar and the p-Ricci tensor as 

E>, . — paia 2 ...a 2 p p \l\ p p.l p P-2 v (A) 

lh {P) J Vi^2...At2 P c ai c a 2 "- C a2 P \ V 

E> . a paaia2...a 2p -i p Hi P V2 p V2 P -i ( c\ 

1> '(P)H >MiM2-M2 P -i °ai c «2 " azp-i > \°J 

Leading to the definitions of the p-Einstein tensor 

G (p)1 = R (p)1 ~ ^ e l r (p) ■ ( 6 ) 

Note that the p-Ricci scalar (|]) coincides with the Euler-Hirzbruch density^ in d = 2p (even) dimen- 
sions, and hence is trivial. In dimensions d < 2p on the other hand, it vanishes by (anti) symmetry 
irrespective if d is even or odd. Thus the upper limit in the summation in (|2|) is P\ = in even, 
and Pi = ==i in odd dimensions d. 

The definition we use for superposed YM hierarchy is 

£ym = £ ^yy r p e Tr F{2pf , (7) 

where F(2p) is the 2p-form p-fold totally antisymmetrised product of the SO(d) YM curvature 2- form 
F{2) 

F{2p) = F^ lfl2 ^ 2p = F[ fllfl2 F^ sfl4 ...F^ 2p _ 1 ^ 2p ] . (8) 

Even though the 2p-form (|5J) is dual to a total divergence, namely the divergence of the corresponding 
Chern-Simons form, the density (^) is never a total divergence since it is the square of one. But the 
2p-form (^) vanishes by (anti) symmetry for d < 2p so that the upper limit in the summation in (|7]) 
is P 2 = | for even d and P 2 = for odd d. 

We define the p-stress tensor pertaining to each term in (0) as 

T« = Tr F(2p), T1T2 ... T2p _ 1 F(2p)/^-^ - ±- 9lw Tr F(2p) nnmn ,F(2p)™~'»> . (9) 

In the present work, we restrict our attention to d — 6 and 8, and in both cases restrict to 
only the two terms p = 1 and 2, both in the graviational and the YM hierarchies, even though in 
these dimensions higher values of P both in (|2]) and (|7|) can be accommodated. This restriction is 
motivated by our desire to study the simplest nontrivial cases. As for our choice of d = 6 and 8, 
these are the only even dimensions in which the system consisting of (0) and (|7|) is nontrivial, and 
our restriction in this preliminary study is motivated by our desire to maximise the analogy with the 
study of Bartnik and McKinnon J^j which is in (even) d — 4. Like in the model of we will choose 
the YM gauge group to be SO(d), with the YM fields taking their values in the chiral representation 
SO±(d) (which for d = 4 reduces to SU±{2)). 

1 This follows immediately from the alternative definition for the p-Ricci scalar in d dimensions: 

£> . _ e v\vi---Vi v H\lil---lia-l v p a\ „02 a d-2 P b 1 b2...b 2p a 1 a2...a,i-2p 0616206364 J}b 2p b2 p 

n {p) ~ fc e Mi e A 1 2--- e A'd-2p t IX V 1 V2 IX V?.Vi--- IX V2 p V2 p ■ 
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3 The classical equations 



In d dimensional spacetime, we restrict to static fields that are spherically symmetric in the d — 1 
spacelike dimensions with the metric Ansatz 

ds 2 = -A(r)dt 2 + Bir^dr 2 + r 2 dtt 2 d _ 2 (10) 

where r is the (d — 2) dimensional spacelike radial coordinate and dQ^-2 is the d — 2 dimensional 
angular volume element. 

We take the static spherically symmetric SO(d) YM field in even d to be in the SO(±)(d) (chiral) 
representation 



The index i = 1,2, ...,d — 1 labels the space dimensions of the static fields, (fi = 1, 2, d) being 
the gamma matrices in d dimensions, sj^ are the generators of the SO(d — 1) subgroup of SO(d), 
the latter being in the chiral representation. We note that (|TTD is not the most general spherically 
symmetric Ansatz for a SO(d) YM connection and is effectively that of a SO(d — 1) field. In this sense 
we could regard our gauge fields to be SO(d — 1) rather than SO(d) fields. In the low dimensional 
case d = 4, the generators S^' reduce to the those of S77+(2) and S77_(2), in which case ( |TT|) is 
indeed general enough. Another difference of d > 4 versus d = 4 is that in the latter case there is a 
duality between the electric and magnetic fields while in our case this duality is absent. Thus we do 
not seek to construct dyonic particle like solutions and have set A = in 

For the system (HD under consideration, namely that with Pi = P 2 = 2, the variational equations 
for the YM and gravitational fields are, respectively 

(e Fn + l -r 2 {F pa , D p (e F^)} = (12) 

^1^(1)^ + ^2^(2)^ = TiT(i)£ + — r 2 T (2) ^ . (13) 

It is easy to systematically extend (0) and (^) to the highest values of P\ and P 2 respectively, in 
any dimension d. 

By virtue of the spherically symmetric YM Ansatz ([TT| ) and the metric Ansatz fllOD , it follows 



that the Gauss Law equation, namely the v = (or v = t) component of (0) is identically satisfied. 
Amongst the v = % components, the v = r also trivialise and the d — 2 equations pertaining to the 
v = 9 a components all coincide to give only one equation for (|TI|). We list explicitly the two terms 
in© 



D,(eF*) = (r^VABw')' - {d - 3)r d ~ 6 J ^(w 2 - l)w (14) 




B 



{F pa ,D,(eF^)} = 3(d-3)(d-A)(w 2 -l)(^(r d - 8 ^(w 2 -l)w')' 



- (d-5)r d - w J^(w 2 -l) 2 w) (15) 



for all a = 1, 2, ...(d - 2). 
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It follows from the metric Ansatz ([H]) that the Einstein tensors defined in ([]) have nonvanishing 
components G^ and G^ e a (a = 1, 2, (d — 2)), which can readily be calculated for any p. 

Here we list explicitly the two cases of interest, p = 1 and p = 2, repectively, 



Gi r 



r (l) 
^ 0. 



2r 

1 

2 



(d-2) 



B' - ~(d-3)(l - B) 
r 

BA' 1 . , 

— --(d-axi-s) 



(16) 
(17) 



iir + V?VA'l-4(^-3)(d-4)(l- J B) 



A 2 \A 
(AB)> 



and 



G (2)t 
Cx(2) r 

G (2)e a 



+ -(tZ-3) , 



_( d _2)(d-3)(d-4)(l- J B) 



B'-_(d-5)(l-B) 



4^3 



(d-2)(d-3)(d-4)(l - B) 



BA' 1 
~X~ 2r 



(d-5)(l - 5) 



(19) 
(20) 



B' 



+ V5) (1 - B )^f 



2r : 



: (d-5)(d-6)(l-B) 5 



(21) 



It follows from the Ansatze C[T0f) and (|TT]) that the stress tensors defined in @ have nonvanishing 
components T(p)*, T( p )£ and Tjj,)^ (a = 1,2, (d — 2)), which can readily be calculated for any p. 
Here we list explicitly the two cases of interest, p — 1 and p = 2, repectively, 



(i)t 

rr-i r 

1 (i)r 

rp da 



2 4 r 2 
1 

2V2 
1 

2V2 



n d (d - 2) 
n d (d - 2) 



2B«/ 2 + (d - 3) 



2Bw' 2 + (d - 3) 



W 2 -l 



2(d-4)Bw' 2 + (d-3)(d-6) 



W 2 -l 



(22) 
(23) 
(24) 



and 



(2)* 



(2)r 



(3!) 2 (d-2)\(w 2 



2 7 r 4 ™ d (d-5)! V r 
(3!) 2 (d-2)! /u> 2 -l' 



(2)e D 



2 7 r 4 (d- 5)! V r 

(3!) 2 (d-3)! f w 2 -l 
2 7 r 4 nd {d-5)\ \~ 



ABw' 2 + (d - 5) 



ABw' 2 + (d - 5) 



w 2 -l 



V 2 -! 



4(d - 8)B«/ 2 + (d - 5)(d - 10) 



w 2 -l 



(25) 
(26) 
• (27) 
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fl24|) and fl2T|) are valid for all a = 1,2, (d — 2). The constants = Tr 2, where the dimensionality 
of the unit matrix is determined by the chiral representations appearing in (|i~lf). 

Substituting (|I^-(|I5|) in ( |I2"D results in one equation, and substituting (|16|)-(f2~T|) in ( |T3"D results in 
three equations. Thus we have 4 ordinary differential equations for the 3 radial functions A(r), B(r) 
and w(r). This necessitates the verification of the consistency of these 4 equations to ensure that the 
system is not overdetermined. We have carried out this consistency check in the case at hand with 
Pi = P 2 = 2 in (Q) in (|?p respectively, for any k,\, k 2 , t±, t 2 and in any dimension d. We found that 
the components (9 a ,9 a ) of (|13"D, 

+ rf (2) £ = nTnll + ^r 2 T {2) e e: . (28) 
are identically satisfied by the two other Einstein equations 

KiG(i)* + k 2 G(2)j = nT(i)* + — r 2 T(2)l (29) 
K\G (i) T r + k 2 G = TiT(i)£ + ^yr 2 T(2)r (30) 

together with the YM equation fll2"|) . Thus the 3 effective constraints on the 3 radial functions are 
the 2 equations (P9"|)-(P0"D plus the YM equation (|T2"D. This check of consistency extends by induction 
to the case of arbitrary Pi , P 2 . 

Finally we express these 3 equations employing the same notation as |]T[ and 0, namely 

A(r) = a 2 (r) N(r) , B(r) = N(r) . (31) 

For the functions w(r), N(r) and cr(r), we have 

2n ((r^crNw')' -{d- 3)r d - 6 a(w 2 - l)w) + 
+ 3r 2 (d - 3)(d - A)(w 2 - 1) ({r d ~ 8 aN(w 2 -l)w')' 



~(d-5)r d - 10 cr(w 2 -l) 2 w)=0 (32) 



Nw' 1 + -(d-3) 

2 \ r 



2 1 \ 2 
ur - 1 \ 



+ 1^(^-3)^-4) 
2 r z 



w 2 -l 



w 2 -l 



Nw' 2 + -(d-5) 

4 \ r 



(33) 



In (^), the function m(r) is defined as 



w 



1 



m {r) = rq l [Kir a ~' 6 {\ - N) + ^ 2 r d " 5 (l - iV) 2 ] 



(34) 



(35) 



which for k 2 = and with d = 4 reduces to the definition of the corresponding function m(r) used 
in ]T|, For k 2 = r 2 == and with (i = 4, equations (|32|)-(|34|) coincide with the ordinary diffrential 
equations of [l], [|. 
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4 Numerical results 



In this section we present the results we obtained for the two cases corresponding to p = 1 and p = 2 
gravity and for the corresponding space-time dimensions d = 6 and d = 8 where solutions of the 
extendend Yang-Mills equations (consisting of the superposition of the p — 1 and p = 2 members 
of the YM hierarchy) exist even in the absence of gravity. (The cases of odd values of d will be 
addressed elsewhere) . When the two mixing parameters T\ , t<i do not vanish we can set them equal 
to particular values without losing generality. This can be done by an appropriate rescaling of the 
overal lagrangian density and of the radial variable. So we take advantage of this freedom and choose 
n = 1,T2 = 1/3 (this simplifies the numerical coefficients in the equations). The limit T\ — > is 
addressed separately. 

4.1 Boundary conditions 

We have solved the above equations with the appropriate boundary conditions for the radial functions 
m(r), air) and w(r) which guarantee the solution to be regular at the origin and to have finite energy. 
These boundary values are 

m(0) 

at the origin, and 

lim a(r) 

r^oo 

at infinity. The condition on a results in Minkowskian metric asymptotically. 
Further analysis of the YM equation (|3^) yields 

»(r - oo) = -1 + ^ (38) 

4.2 p = 1 gravity: Pure Einstein-Hilbert 

We define here a 2 = rid/ (8/«i) so that a 2 = corresponds to the decoupling of gravity (and therefore 
to a flat space). 

4.2.1 d = 6 case 

The flat solution corresponding to the flat case |4j] is characterized by an energy M ~ 7.096 (in the 
conventions used). This solution is smoothly deformed by gravity and exists up to a critical value 
a 2 m « 0.1266. 

When a 2 increases, the mass of the gravitating solution decreases from M ~ 7.096 (flat case) 
down to M m 6.56; the function N(x) develops a local minimum, say N m which becomes deeper 
while gravity becomes stronger and the value <r(0) decreases from one. These quantities are plotted 
in Fig.l and the mass is reported on Fig. 2 

Our numerical analysis strongly indicates that a second branch of regular solutions exists, also 
terminating at a 2 = a^. 

On this second branch the values cr(0) and N m increase monotonically with a 2 . The quantity 
cr(0) reaches the value zero a 2 — >• a 2 , with a 2 ~ 0.07, as illustrated by Fig. 1. 

Accordingly, the metric becomes singular in this limit and no solutions seem to exist on this 
branch for lower values of a 2 . 

The second branch of solutions therefore exist for a c < a 2 << a 2 . The mass of the solution of 
the second branch is larger than the corresponding one on the first branch, as illustrated by Fig. 2. 



w(0) = 1 



lim w(r) 



(36) 
(37) 
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The minimal value of N(r) stays well above zero (N min rs 0.17) so that no horizon is approached. 
The situation contrasts with the case of the gravitating monopole in 3+1 dimensions []5j, ||. In both 
cases (the present and that of || ||) there occur second branches of the solutions after critical values 
the respective parameters (a 2 inourcase) , these bifurcations are quite different in the two cases. The 
gravitating monopole bifurcates into a Reissner-Nordsrtom black hole solution while the solutions 
studied here terminate into a solution presenting a singularity at the origin. Several components of 
the Riemann tensor are indeed polynomials is the field cr(r). 

For scaling reasons, the parameter Ti cannot vanish for finite energy solutions to exist. The 
numerical study of the equation with a varying T\ (with a 2 fixed) reveals that, in the limit t% — > 0, 
the solutions of the first branch converge to the vacuum solution (with mass equal to zero and with 
w(x) = 1) When the solutions on the second branch are examined in the same limit, it appears that 
they stop to exist at a critical value of T\ where the value cx(0) becomes zero. In other words, it also 
terminates into a singular solution for a finite value of T\ and no regular solution seems to exist for 

ri = 0. 

4.2.2 d = 8 case 

The scenario is qualitatively similar to the one of the d = 6 case. Here the flat solution has M ~ 25.6, 
it gets deformed by gravity up to a 2 ~ 0.0022 and the function er(r) is such that <r(0) decreases with 
a 2 . At the same time N(r) develops a local minimum which becomes deeper. Another branch of 
solution exists where both <r(0) and N m i n are lower than their conterparts on the first branch. These 
relevant quantities are plotted in Figs. 3 and 4. On Fig. 5 we present the profiles of the metric 
functions for several values of the parameters. 

The numerical analysis suggests that this second branch persists up to a 2 = and that in the 
limit, lim Q ,2_ +0 cr(0) = 0. As far as our numerical analysis indicates, the value N m tends to a finite 
value in this limit so that there occurs no horizon. 

We further investigated the two branches of solutions in the limit t% — > 0. It appears that the 
first branch solution tends to the vacuum while the second one tends to a non trivial solution with 
a finite mass. 

4.3 p = 2 gravity: Pure Gauss-Bonnet 

This case is not of any particular physical relevance but it is interesting to learn the peculiarities 
of the pure p = 2 gravity especially when the YM term is also the pure p = 2 member of the YM 
hierarchy in d = 8. Such a system supports a self-dual YM solution on double-self-dual gravitational 
background |7J when the metric has Euclidean signature, so it is interesting to contrast with the case 
of Minkowskian signature at hand. 

4.3.1 d = 6 case 

We define the constant /3 4 = 4nd/n2 in this case. Again the flat solution gets deformed by the 
gravitating term and a gravitating solution is formed, characterized by cr(0), N m (see Fig. 6) and by 
the mass (see Fig. 7). 

While both a (0),N m decrease monotonically as functions of /3 2 , the mass first decreases, reaches 
a minimum for (3 2 ~ 0.75 and then starts to increase very rapidly with 1 . 

In the present case, it is very likely that the gravitating solution exists up to a maximal value of 
/3 2 , say (3^ and that, when this limit is approached, N m tends to zero, forming an horizon, while M 
tends to infinity. (Incidentally, tq the value of the radial variable where w(r) crosses zero, exhibits 
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the same behaviour). The numerical integration of the equations for (5 2 > 0.87 becomes difficult and 
we could not evaluate [3 2 with a good accuracy. 

Denoting r m the value of r where the function N(r) attains its minimum, we observe that the 
solutions possess a completely different behaviour in the regions r > r m and r < r m . 

The functions m(r), a(r),w(r) vary essentially on the interval [0,r m ] and attain their asymptotic 
values m(oo), 1, —1 for r ~ r m . They are essentially constant on the interval r e]r m , oof. This feature 
is illustrated in Fig. 8. 

Our belief that the solution stops at some maximal value of ft 2 is reinforced by the fact that our 
attempts to produce a second branch of solutions turned out to be unsuccessful. 

A peculiarity of the equations is that fixing (3 2 and T\ and varying r 2 , is equivalent to a rescaling 
of the radial variable r. Fixing f3 2 and decreasing n — > the numerical results strongly suggest that 
the solution tends to the vacuum. 

4.3.2 d = 8 case 

The pattern of solutions is similar to the ones in the d = 6 case. Here the mass attains its minimum 
at j3 2 ~ 0.55 and again the gravitating solutions seem to exist up to (3 2 ~ 1.2. The relevant quantities 
characterizing the solutions are plotted in Figs. 9 and 10. The behaviour of the solutions for (3 2 
approaching the critical value is illustrated on Fig. 11. 

The absence of a solution for the model consisting purely of the p = 2 members of both gravi- 
tational and YM hierarchies for Minkowskian signature d — 8 at hand, contrasts with the selfdual 
(black hole) solutions in the case of Euclidean signature. (Note that the p = 1 counterpart in 
Minkowskian d = 4 exists, namely the solutions of [1], along with the selfdual solutions.) It would 
be interesting to find out if the situation is different for Minkowskian signature black hole solutions 
in d = 8. 

4.4 p = 1 + p = 2 gravity: Superposed Einstein and Gauss-Bonnet 

We finally considered the equations for the case where the p — 1 and p = 2 gravities are superposed. 
We limited our analysis to the case where the p = 2 gravity coupling is slightly smaller (although 
not too small in order to see the effect of the superposition) than the p = 1 parameter. To fix the 
idea, we choose k 2 = (1/2) aci,. Our conclusion is that for d = 6 and d = 8 the qualitative properties 
of the solutions corresponding to p=l are conserved, namely two branches of solutions are found. 
The maximal value a 2 increases slightly in the mixed case; for instance ~ 0.004 for our choice of 
the parameter in the d = 8 case. 

5 Conclusions 

It appears that the difference of the solutions that we have constructed, from the those given in Jl J 
for the usual EYM system, arises due to the presence of the additional dimensional constants k 2 and 
r 2 . These constants result in the presence of a nontrivial parameter in the equations, for instance the 
one we denoted as a 2 , which contains the gravitational coupling strength. By varying this parameter 
we were able to construct one or several branches of solutions for the cases considered. In this respect 
the system of equations studied is analogous to that of the gravitating monopole where the VEV 
of the Higgs fields also gives rise to an additional dimensional parameter. 

In this prelimary study, we have restricted ourselves to the construction of the simplest static 
spherically symmetric (in d — 1 spacelike dimensions) solutions. Natural extensions of the present 
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work include (i) the construction of solutions with several nodes in the YM function w(r), and 
(ii) the construction of the black hole counterparts of our particle-like solutions. In the latter case 
these pertain to the fully interacting version of the fixed gravitational background solutions discussed 
in @. Another extension, (Hi), concerns the deformation of our solutions due to the inclusion of a 
cosmological constant ||. Finally, (iv) the stability analysis of our solutions can be carried out. 

Acknowledgements This work was supported by an Enterprise-Ireland grant for International 
Collaboration project IC/2001/073. 
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Figure Captions 

Figure 1 The minimal value of the function N and the value <r(0) are plotted as functions of a 2 for p = 1 
gravity in d = 6. 

Figure 2 The mass of the solutions for p = 1 gravity in d = 6. 

Figure 3 Idem Fig. 1 in d = 8. 

Figure 4 Idem Fig. 2 in d = 8. 

Figure 5 The profiles for the metric functions N(r), c(r) are presented for different values of a 2 for the two 
branches for p = 1 gravity in d = 8. The lower branch plot of the a 2 = 0.0008 case is not reported 
because it is very close to the flat solution. Idem [Figure 6] Idem Fig. 1 for p = 2 gravity in d = 6. 

Figure 7 Idem Fig. 2 for p = 2 gravity in d = 6, 

Figure 8 The evolutions of N(r), a(r),w(r), at the approach of the critical value for p = 2 gravity in d = 6. 

Figure 9 Idem Fig. 1 for p = 2 gravity in d = 8. 

Figure 10 Idem Fig. 2 for p = 2 gravity in d = 8. 

Figure 11 Idem Fig. 8 for p = 2 gravity in d = 8. 
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Figure 1 
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